Berry-Phase induced Heat Pumping and its Impact on the Fluctuation Theorem 
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Applying adiabatic, cyclic two parameter modulations we investigate quantum heat transfer across 
an anharmonic molecular junction contacted with two heat baths. We demonstrate that the pumped 
heat typically exhibits a Berry phase effect in providing an additional geometric contribution to heat 
flux. Remarkably, a robust fractional quantized geometric phonon response is identified as well. The 
presence of this geometric phase contribution in turn causes a breakdown of the fluctuation theorem 
of the Gallavotti-Cohen type for quantum heat transfer. This can be restored only if (i) the geometric 
phase contribution vanishes and if (ii) the cyclic protocol preserves the detailed balance symmetry. 
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Understanding and controlling of heat transfer due to 
phonons occurring in low dimensional nanoscale systems 
is both of prime and practical importance Pioneer- 
ing experimental works carried out recently, such as nan- 
otube thermal rectifier nanotube phonon waveguide 
has spawn phononics, i.e. the science and engineering 
of phonons [Ij], as an emerging new scientific discipline 
where heat flow can be manipulated as flexibly as elec- 
tronic current. Although the nonlinear (anharmonic) in- 
teraction has been demonstrated as a crucial component 
0, Q in various functional thermal devices, the heat con- 
trol has heretofore typically been achieved by applying a 
temperature bias, for which in accordance with the sec- 
ond law of thermodynamics - heat flows from "hot" to 
"cold" spontaneously. 

Recent studies show that spontaneous, rare fluctua- 
tions of anomalous heat transfer may occur thus 
being seemingly in apparent violation with the second 
law. Clearly, however, no violation of the second law 
occurs on average. The typical measure of such viola- 
tions is the (small) probability for such anomalous events 
as they emerge from a heat exchange fluctuation theo- 
rem (FT) [6|-|9| . The FT for (nonequilibrium) entropy 
production [Tot and heat flux [7, 8] describes that 
the distribution, Pr{Q), of the heat Q transferred from 
the left (L) bath at temperature Tl to the right (i?) 
bath at Tji over a long time interval r, obeys the rela- 
tion: lim T-Hii[Pr{Q)/Pri-Q)] - Q(/3i?-/3L)/r, where 

r— ^oo 

(3[^j^ = l/ksTL^R- This FT thus shows explicitly that 
heat can transfer spontaneously from "cold" to "hot" 
with finite, although typically with very small probabil- 
ity. In particular, Ref. Q demonstrates this FT in the 
quantum case for heat transfer across a quantum har- 
monic chain coupled with thermal reservoirs. A partic- 
ular challenge that arises is then whether this quantum 
Gallavotti-Cohen type FT remains valid also in the non- 
linear quantum regime beyond the quantum harmonic 
chain limit, and, more generally, whether such a heat- 
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FIG. 1: (Color online) A schematic representation of the an- 
harmonic molecular junction. Quantum heat transfer is gen- 
erated via a dynamics of excitation and relaxation of the local 
single mode. The heat flux J from the center to the right bath 
is deflned as positive. 



flux FT still can be formulated in presence of cyclic time- 
dependent manipulations of certain control parameters. 

In the context of time-dependent manipulations var- 
ious molecular heat pumps have been proposed to effi- 
ciently control heat flux against thermal gradients at the 
nanoscale. In all those cases the system is driven far 
away from equilibrium by use of an external modulation 
imposed on system parameters. For example, a molecu- 
lar model with modulated energy levels, has been found 



to operate as a heat pump [12|. Likewise, a spin sys- 
tem leading to the heat pumping has been studied with 
Ref. 13|. Other schemes investigated pumping of heat 



in electronic nanoscale devices by applying time-periodic 
laser fields [llj. Moreover, Brownian heat motors fueled 
by oscillating temperatures have recently been devised 
as well [l5l . [l6| . Given such time-dependent manipula- 
tions one may therefore scrutinize whether the physics of 
a non-vanishing geometric phase does impact the trans- 
fer of heat under external modulations. If so, what is its 
impact on the existence of a heat-flux FT? 

In this Letter, we shall answer these above mentioned 
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objectives by studying quantum heat transport across an 
anharmonic molecular junction model. We start with a 
system consisting of a molecular junction coupled to two 
thermal baths tl2|, as illustrated in Fig. [TJ The total 
Hamiltonian Htot is composed of the following contribu- 



tions: Htot = Hs 



ttR 



^SB 



g^: system 



Hamiltonian Hs = '^n=o ^■n-\n){n\, with En = nfuvo, 
where we assume that heat transport is dominated by a 
single mode and thus consider a two-level system {N — 2) 
to simulate the strong nonlinearity [l7|. If oo, 
the system reduces to the quantum harmonic case. The 
two thermal baths are represented by sets of indepen- 
dent harmonic modes, i.e., Hg — J^k^kb]. ^,bk.u, with 

ly = L,R, where b\ ^^hk,v are the bosonic creation and 
annihilation operators associated with the phonon mode 
k of bath V. The system-bath interactions is taken to 
be bilinear, i.e., Vgg ~ B^Y^'^Zi ~ 1| + c.c, 

Bv = ')k,v{b\ y -f bk^v)^ where the system-bath inter- 
action is characterized by the phonon spectral function 
V^iuj) = Stt^j. 7^ ^i5(a; — w^.). In the following, we use 
wide-band limit ^^,(10) = r^,. As shown with Ref. 
the limit of fast dephasing and using the Redfield approx- 
imation for weak system-bath coupling, the underlying 
dynamics can be modeled as follows: 

= -piW(fcf_^o + fcf-^o) +PoW(fco''-^i + fco^l^i). (1) 

Here, p„(n = 0, 1) denotes the probability of the molecule 
to occupy the state |n), satisfying po(0 +Pi(*) = 1- The 
activation and relaxation rates read: 

fco%i = r,iv,(wo), fcr^o - r,[7v,(c^o) + 1], (2) 

where N^(ujq) — [e'^"'''^'' — 1]^^ is the Bose-Einstein oc- 
cupation probability. Finally, the steady-state heat flux 
at the right contact (being equal to the heat flux at the 
left contact) is expressed as 



(3) 



where the superscript s means the steady state. The first 
term denotes the energy flux going from the molecule into 
the bath R while the second term provides the opposite 
heat flux from the bath R back into the system. 

Geometric Berry Phase induced Heat Pumping. For 
heat pump operation, the molecular junction connected 
to the two reservoirs is subjected to cyclic parameter 
modulations. This could be realized by imposing a modu- 
lation on either of the following parameters: a;o(i), ^L{t), 
r«(i), TLit), TR(t). Throughout the following, the mod- 
ulations acting on such system parameters are assumed 
to be slow, i.e., we employ adiabatic modulations. Let 
the period of modulation be 7^ = 27r/n. The typi- 
cal frequency for a carbon-carbon bond is 1.4 x 10^'*s~^ 
[l8| . Ti, is around lO^^s"^, according to the measure- 
ment with alkane molecular junction [l9| . The relax- 
ation time for fast thermalization usually is on the order 



of a few fs or ps. Thus, the modulation time scale must 
obey 27r/ri ^ 1 ps. In this way, the assumption of adia- 
batic modulation is valid whenever the driving frequency 
n < 1 THz. 

Of prime interest is the heat flux from the molecule 
into the bath R during the long time span r. This is 
achieved upon introducing the characteristic function for 
the phonon counting field i-e., 2^ 21 1 
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P(0), (4) 
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where P^. (q) is the probability distribution of having heat 
Q — qhujo transferred from the molecule into the bath 
R during time r — >■ 00. Here, 1''' = [1,1], T denotes 
the time-ordering operator, and p(0) = [po(0),pi(0)]^ 
are the intial occupation probabilities. Then, the cu- 
mulant generating function is obtained as: Gix) = 
lim T~^hiZr{x)j which generates the heat current via 

r— >oo 

the relation: J — Hluq dG{x)/d{ix)\^^Q- Denote by 
Ao(XiO the instantaneous eigenvalue of 'H(x,t) with the 
smallest real part and |'0o(X7O) iiVoiXit)]) the corre- 
sponding normalized right (left) eigenvector. The cu- 
mulant generating function takes on the following form, 
being composed of two parts 0,11^, namely: 



Zrix) 



^tQ _ r(gdyn + C/goon 



dtXo{x,t), 



dt{(po\dt\-ipo), 



(6) 
(7) 

(8) 



The flst contribution ^dyn presents the temporal average 
and defines the dynamic heat transfer. This is the only 
term which survives in the static limit. The second, ge- 
ometric part ^goom presents an additional contribution 
caused by the adiabatic cyclic evolution. As we shall see 
it is this part which possesses a nontrivial geometric in- 
terpretation. Let us rewrite ^geom as a line integral over 
the closed contour TZ in the parameter space u: 

t/goom = -7p"^ / du Au, Au ^ {^aiu)\-^\ipo{u)).{9) 

Thus, this is an analog of a Berry phase [2^, which does 
not contain time t explicitly and only depends on the 
geometry of the modulation contour in the parameter 
space u. In the case of two parameters being modulated, 
say ui, U2, using Stokes theorem, we find 



Tp ^ duidu2J^u 1 



(10) 



where S-jz is the integral area enclosed by the contour TZ. 

TuiU2 = (9ui(^o|i9«2'0o) - {du2fo\duiipo) (11) 
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(Color online) (a) The contour map of 
-dJ^TLTRix) / d{ix)\x=o, for Vl = Vr and &Jo = 25 meV. The 
(blue) circle with an arrow denotes the path of two-parameter 
temperature modulations: TL{i) = 200 + 100cos(r2f + 7r/4), 
TR{t) = 200 + 100sin(r2t + 7r/4). The integral area Stz is 
within the circle, (b) Pure Berry-phase induced heat current: 
Jtot = Jgeom {Jdyn = 0). The Straight line is the analytical 
result from Eq. (I14p . while the open circles give the simulation 
results by integrating Eq. ([1} ■ 



is an analog of the gauge invariant Berry curvature [23 1. 

Let us next specify the case that the bath temper- 
atures Ti(t),Tfj(t) are subjected to adiabatic modula- 
tions. Then Eq. ((TT]) yields the Berry curvature in tem- 
perature space, reading: 



^TlTr{x) — -ClCr 



2ism{x)TLTR{rL+TR) 



(V^2 + 4£))3 



(12) 



where kBl3lhujoe'^"'"^'>N^, K = Tl{1 

Tr{1 + 2Nr), D = TLVRNLNR{eP^'''-«e+^ + ef'^^'^^e^^) 
with e±x = e^*-^ — 1. Upon substituting this Berry cur- 
vature into Eq. ([T0| . the total heat flux emerges as: 



9[^dyn(x) +^geom(x)] 



Jd 



yii \ "^georn; 



^0 ,SlTr{Nl-Nr) 



J gr, 



where 
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d{ix) 



2ClCrTlTr{Tl+Ti 



(13) 
(14) 

(15) 



The dynamic part Jdyn just coincides with the tempo- 
ral average of the heat flux obtained from J = J{t) in 
Eq. dS]). The geometric part Jgoom is the additional heat 
flux that results from the nontrivial Berry phase effect. 
The ratio of this geometric heat flux and the dynamic one 
is typically about Vl/Ti,. To avoid that Jgeom is masked 
by Jdyn, we choose a symmetric molecular junction with 
= Tr, and modulate TL{t),TR(t) as indicated by the 
circle contour in Fig. [5^. Then one finds that Jdyn = 



and Jgcom 7^ 0, see Fig. [2Jd, so that the Berry phase in- 
duced Jgeom dominates the heat transport. This is the 
case for which the geometric phase effect on heat trans- 
port is distinctly experimentally detectable. As a main 
finding we have that the Berry phase effect acts as a heat 
pump, providing an additional heat flux across the molec- 
ular junction even though on average no thermal bias acts 
and the system is symmetric. Note also, distinct from the 
irreversible heat flux Jdyn, Jgoom is time- reversible, i.e., 
under the time-reversed modulation (t — > —t) the Berry- 
phase induced heat fiux just reverses sign. 

Fractional Quantization of Phonon Response. Re- 
markably, we find a fractional quantized phonon response 
for large temperature driving: the integral in Eq. (jl4p can 
be rewritten asj^ dNLdNR2TLrR{TL + Tr)/K^ = 
1/4, yielding [2a |: 



Jgcom — -^fvjJo/Tp. 



(16) 



This 1 /4 fractional quantized geometric phonon response 
is robust since it does not depend on the specific values 
of Hloq^T i^,T R. It means that the geometric phase effect 
caused by two bath temperature modulations is able to 
pump maximally on average one phonon ?iajo per four 
cycles. 

Impact on Heat-Flux Fluctuation Theorem. Besides 
the dynamic t^dyn, ^gcom will generally not only con- 
tribute additionally to the average heat transfer but also 
impact the higher moments of the heat current (such 
as the phonon counting statistics) and other heat trans- 
port characteristics as well. In the following, we study 
its impact on the heat-fiux fluctuation theorem. Before 
doing so, let us address first the static situation with 
^/gcom(x) = 0, yielding 



Q{x) = ^dyn(x) = -Ao(x) 



-K + ^/lO+AD 



(17) 



Then, Aq obeys the Gallavotti- Cohen (GC) symmetry 
di] (and alike for tJdyn(x), ^^(x) and 2t(x)), reading 



Ao(x) = Ao(-x + */?*), 



(18) 



where /3* = Hiko^ikf-^o) I {ko-^iki^o)]- In virtue of 
Eq. ([2|), yielding the detailed balance relation fcg^^ — 
k'^^oe'l^"'"^°, we find that = hujoil3R - Pl)- Via an 
inverse Fourier transform of Eq. (j?]), this GC symmetry 
results in the quantum FT of heat transport for an an- 
harmonic molecular junctions, reading with Q = qHuiq: 



lim ihr[4M_] 



= QWb. - Pl)/t, 



(19) 



which precisely coincides (without any correction) with 
the result for the quantum harmonic chain Q. This FT 
gives the probability of observing spontaneous "second 
law violation": Assume Tl < Tr, i.e. (3* < 0; the upper 
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bound to observe the "violation" for spontaneous heat 
transfer from (left) cool to (right) hot is estimated as: 
dqPriq) = dqPr{~q)e'f^' < e'=^\lt indicates that 
in absence of external modulations, the probability of at 
least c phonons (or net energy chuo) transporting against 
the thermal bias is nonvanishing detectable, although de- 
caying exponentially. 

For the time-modulated system the GC symmetry 
ceases to hold when tjgcom(x) 0. For example, in the 
case of cychc temperature modulations TL{t) and Tii{t), 
the Berry curvature ^TlTr{x) contains the factor sin(x), 
which explicitly breaks the GC symmetry of 5goom(x)j 
and alike for g{x) and Z^(x)- Thus, the FT Eq. ^ 
becomes violated as a consequence of a geometric phase 
induced breakdown of GC symmetry. Moreover, even 
for parameter modulations yielding Qgcomix) — Oi and 
with time-dependent (3* — ^ fi*{t)^ the GC symmetry for 
^'dyn(x) = -Tp~^ Jq" dtXo{x,t) generally cannot be re- 
covered, despite Xo{x,t) — Ao(— x + */3*(^), i)- 

Interestingly, we find that for time modulations of 
the system-bath couplings ri^{t),r]i{t), the detailed bal- 
ance relation fco^i/fci-^o ~ e"^"'*'^'' remains intact, thus 
providing a vanishing Berry curvature -^ri,rjj(x) = 0- 
Meanwhile, with the resulting time-independent f3*{t) = 
/?*, one finds that the GC symmetry of tJdyn(x) = 
~Tp~^ /o^" dtXo{Xjt) still holds. Consequently, we obtain 
a vanishing Berry-phase induced heat pumping and, sur- 
prisingly as well, also no violation of the FT, no matter 
how TL{t) and rfi{t) are modulated. 

In summary, through investigating heat transport 
across an anharmonic molecular junction by applying 
cyclic two-parameter modulations, we find that the sys- 
tem generally undergoes, apart from dynamic pumping, 
also a Berry phase induced heat pumping. This geo- 
metric contribution exhibits a robust fractional quan- 
tized phonon response. Furthermore, the quantum FT 
for heat transport in presence of a static temperature bias 
holds true in the anharmonic case as well. The presence 
of the geometric phase, however, violates the heat-flux 
FT. Only in situations of vanishing Berry curvature and 
restoration of detailed balance symmetry can the validity 
of the FT be recovered. 

Although our present work did focus on the adiabatic 
regime, it likely can be extended to the case of a non- 
adiabatic geometric phase [25], and maybe also for non- 
cyclic modulation schemes in the spirit of ^22!]. Because 
the geometric phase has profound effects on material 
properties 23| we hope that our present findings do in- 
vigorate others to undertake related studies aimed at un- 
covering intriguing novel geometric phase induced ther- 
mal effects (such as thermoelectricity) which will enrich 
further the discipline of phononics. 
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SUPPLEMENTARY INFORMATION OF "BERRY-PHASE INDUCED HEAT PUMPING AND ITS 

IMPACT ON THE FLUCTUATION THEOREM" 



The Generating Function with the Phonon Counting Field and associated Geometric Phase 

In this section, we derive the geometric phase effect from time evolution of the generating function with phonon 
counting fields in more detailed form. For our system, Pt (q) denotes the probability that within time t the net number 
of phonons transferred from the molecule into the bath R is q. We then split this probability into two parts, namely, 
Ptiq) = Pt{0,q) + Pt{l,q), where Pt{0,q) {Pt{l,q)) denotes the probability that having q net phonons transferred 
from the molecule into the bath R, within time t, while the molecule is dwelling on the low "0" (high "1") energy 
level at time t. The time evolution of our system can be described as: 

^Pt(0,<z) = -(fco^i+A:o^i)Pt(0,<z) + fcf_,oPt(l,q) + /cf:,oPt(l,g-l), 

j^Pt{l,q) = fco^-^iPt(0,g) + fco''-^iPt(0,(Z + l)-(fcfl^o + fcf_^o)Pt(l,(7). 

By multiplying the factor e^"^^ on both sides of the equations and summing over q, we obtain the time evolution 
equation of these individual generating functions: 



(SI) 



where 



l,R 



l*(x,t)> 



Therefore, the characteristic function Zr{x) emerges as: 

■DO 



Jo 'H{x,t)dt 



l*(x,o)). 



(S2) 



q— — oo 



At initial t — 0, the net number of transferred phonons is 0, thus |^'(x, 0)) = [poit — Q),pi{t — 0)]"^. Note that the 
time evolution equation Eq. (SI) for generating functions including the counting field exists for systems described by 
such master equations generally [1]. 

In the following, we are going to derive the geometric phase from the time evolution of the generating function (we 
omit the counting field variable x for the sake of better readability). We stress that the time-dependent operator H{t) 
is not a genuine, hermitian Hamilton operator. However, we still can decompose it as 

Hit)\ibnit)) = A„(t)|?/;„(i)), {^„{t)\Hit) = ((p„(t)|A„(t), 

where A„(i) denotes the instantaneous eigenvalue of 'H(i) and \ipn(t)) ((</?n(i)|) the corresponding normalized right 
(left) eigenvector, satisfying {ipmit)\4'nit)) — Smn- Similar to Berry's reasoning Q, we expand |\E'(i)) as 



|*W> =X^a„(i)exp(^-^ Xn{t')dt'^ \Mt))- 
By substituting Eq. (S3) into the time evolution equation Eq. (SI), we find 

^d„(i)exp (^-^ Xn{t')dt'^ IMt)) =X!""W^^p("/ K{t')dt'^ \Mt))- 
Upon the multiplication of the left eigenvector ((/?m(t)| and observing that {ipm{t)\^n{t)) = < 



(S3) 



(S4) 



we obtain 



amit) = -a,nit){ip„^it)\^P„^it)) - ^ a„(i)exp ( - / [X„{t') - \rnit')]dt' ) {^„,{t)\i:nit)). (S5) 
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Note An are complex- valued eigenvalues. The long time behavior of the system is governed by the eigenmode to = 0, 
whose eigenvalue Aq possesses the smallest real part. Consequently, by neglecting the last term within our adiabatic 
approximation, one obtains: 



«o(t") = exp 

Taking into account the adiabatic cyclic evolution over a long-time period 7J,, we end up with 

= lt|*(T)) « exp j^' dt [Ao(t) + (^ol^o)] j ao(0)lt|^o(T)), 



(S6) 



(S7) 



which tells us that the cumulant generating function Q (%) = lini t ^InZr (x) = Gdvn + Gseom contains two contri- 

r— ^oo 

butions, one originating from the dynamic phase factor and the other from the geometric phase factor (the boundary 
value contribution: In [ao(0)l'''|^/'o(T))] becomes negligible in the long time limit): 



Gdyn = -Tp ^ / diAo(x,t), 
Jo 



(S8) 



/ dt{^o\dt\i^o) 



(S9) 



These two expressions of Gdyn and t?gcom coincide precisely with those obtained in Ref . [3] , although derived therein 
using a different approach. 



Fractional Quantized Phonon Response 

Now, we give a detailed explanation of the physical picture of the observed 1/4 fractional quantized phonon response, 
in particular, the cause of the Berry phase effect induced by temperature modulations of the two heat baths to transfer 
on average 1/4 phonon hwo per one modulation cycle. 

Mathematically, when the integral area S-jz due to the two parameter modulation encloses a large part of the 
area around the maximum of — 9JTiTjj(x)/t^(*x)lx=0: O'"! Put differently for strong temperature driving temperature 
comprising near zero temperature up to large values beyond fkoo, as sketched by the dashed path in Fig. [3^, the 
integral in Jgcom can be recast as 

^gcom = — J J dNLdNu — = -hwo/Tp, (SIO) 

The underlying physical mechanism is as the following: As shown in Fig. [3^, the trajectory of this so obtained 
temperature modulation follows the path: (1) (2) (3) (4) (1), where (1) (Tl = 0,Tr = 0), (2) (Tl = 
0,Tr -> oo), (3) (Tl oo,Tr oo), (4) {Tl oo,Tii ^ 0) as detailed in Fig. For the two-level system under 
consideration, only for the parameter set around (1) does the system fully occupy the lower level; at the other three 
sets (2), (3), (4), the system occupies the upper level and lower one with equal probability. The level transitions 
during the course of the temperature modulation are illustrated with Fig. [3J;. Usually, the transition from the lower 
(upper) level to the upper (lower) contain two parts of contributions: one is from the left bath and the other is from 
the right bath, respectively (see Fig. 1 in the text). However, for the transition from (1) to (2), the temperature of 
bath L is increased from to oo so that only the up transition from bath L is excited, which is indicated as the red 
arrow; while for the transition from (3) to (4), Tl is modulated near zero and Tr keeps extremely high so that only 
the transitions from the R bath become allowed. Those are indicated as blue arrows (see Fig. |3t). Besides, black 
arrows indicate that the system stays intact and gray arrows depict the transitions occurring at either L or R bath 
with equal probability. Note that only the up transition at L followed by one down transition at R counts for the 
positive energy transport from L to R. 

The transitions during temperature modulations can be further decomposed into eight individual paths, as illus- 
trated in Fig. U Therefore, we can count the amount of energy transport for each path individually. For example, in 
path (a), the up transition is excited at L during (1) — > (2), and then the system keeps in the upper level until the 
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FIG. 3: (Color online) A schematic representation of the temperature modulation cycle and corresponding level transitions. 
(a)The contour map of —dJ-Ti^Tnix) / '^i^x)\x=Oj for = Fh and hwo = 25 meV. The (red) dash line with an arrow denotes 
the trajectory of temperature modulations. (b)The temperature list of state (1)(2)(3)(4). (c) The corresponding transitions 
during cycle (1) — >■ (2) — ^ (3) — >■ (4) — >■ (1). The red arrow denotes only the up transition from bath L is excited; while the 
transitions only allowed at the R bath are indicated by blue arrows. Besides, black arrows indicate that the system stays intact 
and gray arrows depict the transitions occurring at both L and R bath with equal probability. 



down transition happens during (4) — ^ (1). Since the down transition exhibit a splitting, occurring at L or R, this 
path contributes 1/2 phonon for the energy transport. In path (b), the up transition at L is first excited and then 
during (3) — > (4) the down transition occurs at R, so that one phonon is transported during the whole modulation 
cycle. The amount of energy transport in other paths can be accounted likewise. 

Finally, the total amount of phonon transport during one temperature modulation cycle is expressed as the average 
of contributions from these eight routes: 

1/2 + 1 + + 1/2 + + (-1/2) + + 1/2 _ 1 

8 " 4- ^^^^> 

In other words, after one complete modulation cycle, the system returns to its "original" state but with hajo/A energy 
given away. This effect has a geometric interpretation and can be utilized to act as a heat pump. 

In summary, this scenario makes it plausible that the Berry phase effect caused by two-parameter temperature 
modulations is able to induce on average a heat transfer across the molecular junction in the amount of 1/4 phonon 
hajQ per modulation cycle. This 1/4 fractional quantized phonon response is robust in the sense that it does not 
depend on the specific values of huo^TL^Tn. 
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FIG. 4: (Color online) Decomposition of the level transitions into eight individual transition paths. The amount of phonon 
transport contributed by each path is indicated. 



